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The multi-scale entanglement renormalization ansatz (MERA) can be used, in its scale invariant 
version, to describe the ground state of a lattice system at a quantum critical point. From the scale 
invariant MERA one can determine the local scaling operators of the model. Here we show that, 
in the presence of a global symmetry Q, it is also possible to determine a class of non-local scaling 
operators. Each operator consists, for a given group element g £ Q, of a, semi-infinite string Fg with 
a local operator (fi attached to its open end. In the case of the quantum Ising model, Q = 1ii, they 
correspond to the disorder operator /i, the fermionic operators and and all their descendants. 
Together with the local scaling operators identity I, spin a and energy e, the fermionic and disorder 
scaling operators 7/), 1/) and /i are the complete list of primary fields of the Ising CFT. Thefore the 
scale invariant MERA allows us to characterize all the conformal towers of this CFT. 



PACS numbers: 03.67.-a, 05.50.+q, 11.25.Hf 

The multi-scale entanglement renormalization ansatz 
(MERA) [U [2] is a tensor network introduced to effi- 
ciently represent ground states and low energy subspaces 
of quantum many-body systems on a lattice. It is based 
on a real-space renormalization group (RG) technique 
known as entanglement renormalization [I] , that employs 
unitary tensors (disentanglers) to remove short-range en- 
tanglement from the system at each RG iteration. The 
removal of entanglement is a key difference with other 
real-space RG techniques, such as Wilson's ground break- 
ing numerical RG (NRG) for the Kondo problem [3j or 
White's extremely succesful density matrix RG (DMRG) 
0] for arbitrary one-dimensional systems. At a fixed 
point of the RG flow, it produces a representation that 
is explicitly scale invariant: the scale invariant MERA. 
This ansatz is characterized by only a small number of 
tensors and can be used to describe systems with topo- 
logical order (at a non-critical RG fixed point) [S^; as well 
as continuous quantum phase transitions (corresponding 
to a critical RG fixed point) [H El EHTO]. 

Here we shall be concerned with the characterization of 
a (scale invariant) quantum critical point with the scale 
invariant MERA [Tl [SHTU]. Evidence that the scale- 
invariant MERA is capable of describing critical ground 
states was first presented in Ref. [1] for the quantum 
Ising model and in Ref. [5] for non-interacting systems 
of fermions and bosons. On the other hand it was argued 
that this ansatz naturally reproduces two important as- 
pects of critical ground states: the logarithmic scaling for 
the entanglement entropy of a block of contiguous sites 
(in one-dimensional systems) and the power-law decay of 
correlations 0. The latter was seen to follow from the 
fact that a two-point correlator C (31,82) between two 
points separated a distance r = |si — S2I is obtained af- 
ter 0(log(r)) applications of a fixed superoperator that 
introduces a constant factor z < 1 after each application, 
and therefore 

C2isi,S2)^z'°Sir) =r-'^^ q=-\0gZ. (1) 



This result was formalized in Ref. [7] by relating the 
possible values of the factor z with the eigenvalues of 
that superoperator, and by identifying the scaling opera- 
tors of the theory with the corresponding eigenoperators. 
A connection between the scale invariant MERA and the 
conformal field theory (CFT) underlying a quantum crit- 
ical point was then established in Ref. [5j, including a 
way to extract the conformal data: central charge, pri- 
mary fields, and their scaling dimensions and operator 
product expansion (OPE) [12,. However, the analysis of 
Refs. [71 IH] was only concerned with local scaling opera- 
tors. For instance, for the quantum Ising model, Ref. [5] 
identified the primary fields identity I, energy e and spin 
a, as well as some of their descendants, all of which were 
expressed as an operator acting on two contiguous sites 
of a coarse-grained lattice. Instead, non-local scaling op- 
erators were not considered. One reason is that entan- 
glement renormalization, being based on locally coarse- 
graining the system, has a computational cost that grows 
exponentially with the size of the support of the opera- 
tors under consideration. 

In this paper we show that the scale invariant MERA 
can be used to characterize a whole class of non-local 
scaling operators of a critical lattice model. We consider 
a quantum spin chain whose Hamiltonian H is invariant 
under a symmetry group G, 

TgHTl= H, yg e g, (2) 

where Tg = ■ ■ ■ Vg (E) Vg Vg ■ ■ ■ is an infinite string of 
copies of a matrix Vg, with Vg a unitary representation 
of g. We shall see that, by incorporating the symme- 
try Q into the MERA, it is possible to study non-local 
operators that have a semi-infinite string of Vg's. Impor- 
tant examples of such non-local operators are the disorder 
operator fi and the fermionic operators ^ and ip of the 
quantum Ising model, which are associated to the low 
energy spectrum of a chain with anti-periodic boundary 
conditions. 
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FIG. 1. (i) Scale invariant MERA, characterized by a dis- 
entangler u and an isometry w that are copied throughout 
the ansatz. (ii) These tensors are isometric, meaning that 
u^u = /, lo^w = I. (iii) We choose u and w to be symmetric, 
Eq. [sj (iv) As a result, the infinite string Eg commutes with 
a layer of disentanglers and isometries. In other words. Eg is 
invariant under coarse-graining. 




FIG. 2. (i) Coarse-graining of a non-local operator o"^ = Eg (g) 
o. (ii) Most of the string ES of Vg's commutes with the coarse- 
graining thanks to Eq. ^ (iii) Then we can remove most 
of disentanglers and isometries using Fig. [ijfi). (iv) o' is 
defined in term of o, u, w and Vg. (v) Scaling superoperator 
Sg, o' — Sg{o), for the local part of a non-local operator Og 
in Eq. [6] see Eqs. [7]|8] Notice the average over the three 
possible ways in which o can be coarse-grained. In the case 
of (7 = I, we have Vg = I, so that that Oj is simply a local 
operator and we recover the 'usual' scaling superoperator 5 
of Eqs. rais] (see Fig. 1 of Ref. [8] for further details). 



Local scaling operators. — Recall that the scale invari- 
ant MERA is made of copies of a unique pair of bulk 
tensors, namely a disentangler u and an isometry w, dis- 
tributed in layers according to Fig. [iji). In the presence 
of the symmetry (pi), we choose these tensors to be in- 
variant under Q [Tlj . 

{Vg ^ Vg) U {Vg ® Vg)^ = U, 

{Vg ^Vg^ Vg) W {Vg)^ = W (3) 

where Vg acting on different indices may actually denote 
different (in general, reducible) representations otQ. The 
layers of disentanglers and isometries define a real space 
RG transformation and a sequence of increasingly coarse- 
grained lattices {L, L', L", • • • }■ Under coarse-graining, 
a local operator o transforms according to the scaling 
super-operator S of Fig. [2]jv) for 5 = 1, 

— > — > o ■ ■ ■ (4) 

The scaling operators (j)a and scaling dimensions are 
obtained from the eigenvalue decomposition of the scaling 
superoperator 5 [II [8] , 



S{ 



A. 



(5) 



Non-local scaling operators. — In this work we consider 
the coarse-graining of non-local operators of the form 



■Va®V„®Vn 



(6) 



where T^ is a semi-infinite string made of copies of Vg 
and o is a local operator attached to the open end of F^. 
Notice that, under coarse-graining, is mapped into 
another non-local operator o^' of the same type. 



F>o 



(7) 



since the semi-infinite string Tg commutes with the 
coarse-graining everywhere except at its open end, as il- 
lustrated in Fig. [2] where we exploit that the disentangler 
u and isometry w have been chosen to be symmetric, Eq. 
[3] In other words, we can study the sequence of coarse- 



grained non-local operators 



by just 



coarse-graining the operator o with the modified scaling 
superoperator Sg of Fig. [2j 

— > a — > ■ ■ ■ (8) 
In particular, by diagonalizing this scaling superoperator, 

Sg{4'g^a) = ^g,a<t)g,a , ^g,a = ^ loga Ag,Q , (9) 

we obtain non-local scaling operators i/)^ „ of the form 



(10) 



Notice that for 5 = I we recover the local scaling opera- 
tors 0a of Refs. [71 [5] . 
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Quantum Ising model. — As a first example, we use the 
above formalism to identify the non-local operator con- 
tent of the Ising CFT starting from the Ising quantum 
spin chain, as described by the Hamiltonian 

oo 

i/,..„,^ J2 {X{r}X{r + l) + Z{r + l)), (11) 

r— — oo 

where X and Z are Pauli matrices. This model preserves 
parity, G = ^2, so that g E {+1, —1}, with V+i = I and 
y_i = Z, and 
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TABLE I. Scaling dimensions of a few non-local operators of 
the quantum Ising model. The conformal towers of -0 and tp 
have identical scaling dimensions. 



T-i = 



Z. (12) 



m— — OQ 



Each index i of tensors u and w decomposes as i = 
(p, ttp), where p labels the parity (p = 1 for even par- 
ity and p = —1 for odd parity) and ap labels the dis- 
tinct values of i with parity p. Then the tensors u, w 



are chosen to be parity preserving, e.g. 



if 



p{ii)p{i2)piji)pij2) = —1- An operator O acting on the 
spin chain has parity p if (F-i) O (r_i)^ = p O. 

We have used the algorithm of Refs. [H [13] to obtain 
a scale invariant MERA approximation for the ground 
state of i?i3i„g, with a computational effort that scales as 
0{x^X^) with the dimension x (and x) of the lower (and 
upper) indices of the disentangler u. The present results 
correspond to x = 36 and x — 2Q and required one week 
on a 3 GHz dual core desktop with 8 Gb of RAM. The 
scaling superoperators Si and S-i were diagonalized in 
each parity sector. The resulting non-local scaling oper- 
ators are of the form 



(13) 



Table [T] contains a few scaling dimensions extracted from 
5-1. The second and fifth columns are for scaling op- 
erators with even and odd parity, respectively, and re- 
produce the exact results with several digits of accuracy. 
Fig. [3] shows scaling dimensions for both local and non- 
local operators. Local scaling operators with even parity 
form the two conformal towers |12) of the primary fields 
identity I and energy e of the Ising CFT, whereas those 
with odd parity form the conformal tower of the primary 
field spin a. Non-local scaling operators with even par- 
ity form the conformal tower of the disorder operator /z, 
and those with odd parity are organized according to two 
towers corresponding to the fermion operators ■0 and tp. 

We have also computed the coefficients Cap-t of the 
operator product expansion (OPE) |12) for all primary 
fields, by analysing three-point correlators as explained 
in Ref. p.. Notice that a three-point correlator 
('f'giai4'g2p4>g3'y} ^111 vaulsh unlcss (i) the product of par- 
ities of the three operators is +1 (since the ground state 
is invariant under parity) and (ii) 315253 = I G 5 (since 
otherwise the product 4'giai'l>g2p'l^g3-y ^ non-local op- 
erator o'^, which must decompose as a sum of non-local 
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FIG. 3. A few scaling dimensions of local (left) and non- 
local (right) scaling operators of the quantum Ising model, 
organized in its six conformal towers. 



scaling operators 0^, and {(p'') — since all non-local 
scaling dimensions A_i^q are larger than zero, so that 
(o^) — 0). Table [lT| shows a numerical estimate of all 
non- vanishing OPE coefficients Ca^-y- Again, the results 
match the exact solution with several digits of accuracy. 
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TABLE II. OPE coefficients for the local and non-local pri- 
mary fields of the Ising CFT. 

Thus, not only have we been able to identify the en- 
tire field content {I, e, cr, t/j, -0, /j,} of the Ising CFT from a 
simple and rather unexpensive analysis of a quantum spin 
chain, but we can now also identify all possible subsets 
of primary fields that close a subalgebra by inspecting 
Tabic |TTj Indeed, it follows that we have the following 
fusion rules 

exe = I, o-X(T = I + e, CTXe = CT, (14) 
/iX/i = I + e, fi X e — fi, (15) 
0xV = I, ■0x-0 = I, (16) 
?/)Xi/j = e, ■0X6 = 0, 05xe = -0, (17) 
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(as well as other, such as cr x /i = + ■0, etc) from where 
we see that {I, e} and {I,e,cr} close subalgebras of local 
primary fields, whereas {I, e, fj,} and {I, e, tp, ■tp} close sub- 
algebras that contain both local and non-local primary 
fields, where locality is relative to the spin variables. 
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FIG. 4. (i) Some scaling dimensions for local operators of 
the quantum XX model. [Sectors with particle numbers +\n\ 
and — |n| yield the same scaling dimensions.] (ii) Some scaling 
dimensions for non-local operators with Vg — Z. (iii) Scaling 
dimensions Ag^a as a function of 9, see Eq. |19| The scaling 
dimensions for n = 0, 1, 2 appear to only differ by a shift. 



The optimization of a scale invariant MERA with x = 
54 and x = 32 took one week (by exploiting the block 
structure of the tensors [H]). For several values of 6* e 
[0,27r), we diagonalized the scaling superoperator Sg in 
the particle number sectors n = 0, ±1, ±2, • • • . Fig. Qi)- 
(ii) show the resulting scaling dimensions Ag^a for = 
and TF — that is for Vq = 1 (local operators) and = Z, 
which also appear clearly organized in conformal towers. 
Finally, Fig. Qiii) shows the scaling dimensions as a 
function of 6. They are seen to accurately approximate 
the expression (denoted 'exact' in Fig. [dj^iii)) 



Ao., 



2tt 



q = 0,±l, (20) 



which is consistent with previous results ^Mj. Up to a 
shift and a rescaling factor, the scaling dimensions Ag^a 
reproduce the low energy spectrum of the XX chain with 
twisted boundary conditions with twisting angle 9. 

In summary, we have explained how to use the scale in- 
variant MERA to characterize non-local scaling operators 
of a critical quantum spin chain. For the quantum Ising 
model, we have identified all non-local primary fields and 
obtained remarkably accurate estimates of their scaling 
dimensions and OPE coefiicients. For the quantum XX 
model, we have obtained continuous families of non-local 
operators associated to twisted boundary conditions. 

Support from the Australian Research Council (APA, 
FF0668731, DP0878830) is acknowledged. 



Quantum XX model. — As a second example we study 
the quantum spin chain with Hamiltonian 

oo 

i7xx= {X{r)X{r + l)+Y{r)Y{r + l)), (18) 

r— — oo 

where X and Y are Pauli matrices. This model is invari- 
ant under rotations Vg — e^*^^/^ on all spins. Therefore 
Q = U{1), group elements g can be labeled by an angle 
6 e [0,27r), and 

OO 

Tg i/xx rl = i/xx, Fe = (g) Vg. (19) 

m— — oo 

To simplify the analysis, we regard each site as containing 
two spins, so that the on-site z-component of the spin can 
take the values and ±1. Then each index j of a tensor 
decomposes as i — (n,a„), where n S Z is the 'particle 
number' (z spin component) and a„ labels the distinct 
values of i with particle number n. Tensors u and w 
are chosen to be invariant under U{1), e.g. uj^J^ = if 
n{ii)+n{i2) n{ji)+n{j2). An operator O acting on the 
spin chain has particle number n if {rg)0{Tgy = e~™^. 
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